Transport properties of anyons in random topological environments 
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The quasi one-dimensional transport of Abelian and non-Abelian anyons is studied in the presence 
of a random topological background. In particular, we consider the quantum walk of an anyon that 
braids around islands of randomly filled static anyons of the same type. Two distinct behaviours 
are identified. We analytically demonstrate that all types of Abelian anyons localise purely due to 
the statistical phases induced by their random anyonic environment. In contrast, we numerically 
show that non-Abelian Ising anyons do not localise. This is due to their entanglement with the 
anyonic environment that effectively induces dephasing. Our study demonstrates that localisation 
properties strongly depend on non-local topological interactions and it provides a clear distinction 
in the transport properties of Abelian and non-Abelian statistics. 



PACS numbers: 05.30.Pr, 03.65.Vf, 05.40. Fb 

In systems with physics constrained to two dimen- 
sions, point like particles named anyons can occur which 
have more general statistics than bosons or fermions [1]. 
Beyond mere possible existence they were found to be 
a good description for low lying quasi-particle excita- 
tions of fractional quantum Hall systems, Majorana edge 
modes of nanowires, and they exactly describe excita- 
tions in various strongly correlated two dimensional spin 
lattice models [2]. Recently there has been experimental 
progress in preparation and control of systems capable of 
exhibiting topological order with the goal to observe any- 
onic statistics [3, 4]. This is further motivated by the dis- 
covery that braiding some types of non-Abelian anyons 
provides for naturally fault tolerant quantum computing 
[2, 5]. As we are not yet able to manipulate anyons in- 
dividually, it would be beneficial to reveal their exotic 
braiding properties on macroscopic scale. In particular, 
we are interested in the transport properties of anyons 
and their possible localisation, that can have direct ob- 
servable consequences. 

Transport properties of anyons in ordered backgrounds 
were studied in [6]. There a discrete-time quantum walk 
was considered of one mobile anyon braiding around a 
canonically ordered set of static anyons positioned on 
a line. While for Abelian anyons the dispersion of the 
walker is quadratic, as in usual quantum walk [7], the 
non-Abelian anyons have richer behaviour induced by 
their braiding properties. Specifically, it was shown in 
[ ] that the Ising non-Abelian anyons exhibit asymptoti- 
cally a linear dispersion relation, characteristic to classi- 
cal random walks. Entanglement resulting from braiding 
the mobile non-Abelian anyon around the static ones is 
sufficient to suppress quantum correlations responsible 
for the quadratic speed-up. 

In this letter, we investigate the role of disorder on the 




FIG. 1: The quasi one-dimensional quantum walk of an anyon 
braiding counterclockwise around islands filled with a random 
number of static anyons of the same type. The islands, de- 
noted with dashed circles, are canonically arranged on the 
line. The possible positions of the walker are denoted by dot- 
ted circles placed in between the islands. 



propagation of both, Abelian and non-Abelian, anyons. 
It has been known for more than five decades that ran- 
domised local potentials can suppress diffusion of quan- 
tum particles — a phenomenon known as Anderson lo- 
calisation [9]. This mechanism is based on randomisa- 
tion of phases that correspond to individual particle his- 
tories and consequent destructive interference. Here we 
consider a discrete-time anyonic quantum walk in disor- 
dered topological background. In particular, we consider 
a walker that braids around islands canonically arranged 
on a line, where the number of static anyons at a given 
island assumes a random value. For the Abelian anyons 
this causes the walker acquiring random discrete phases 
which, as we demonstrate, leads to localisation. For 
the non-Abelian ones we resort to numerical treatment. 
A continuous-time analog, studied numerically with the 
matrix product state approach, concludes that Ising non- 
Abelian anyons do not localise. These generic behaviours 
provide a very clear distinction between Abelian and non- 
Abelian anyons. 

The discrete-time model:- Our setup consists of n "is- 
lands" canonically ordered on the surface and labeled by 
index s, as shown in Fig. 1. The s-th island is occupied 
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by m s static anyons (m s > 0) and the configuration, 
represented by vector m = (mi, . . . , m n ), is supposed 
to be fixed during the course of the walk. Anyons are 
labeled within an island from left to right by an index 
i s = 1 , . . . , m 8 . The mobile walker anyon hops between 
neighbouring sites winding counterclockwise around the 
islands. Hence our system is quasi one-dimensional. We 
denote possible walker's spatial positions also by s, with 
the convention that position s lies between islands s — 1 
and s, as shown in Fig. 1. The distance between sites 
is set to unity, though for purely topologically interac- 
tions the distance scale is irrelevant. Hopping direction 
is controlled by the coin state: |0) moves the walker to 
the left, |1) to the right. The total Hilbert space decom- 
poses as H = "Hspaco <8> "Hcoin ® 'Hfusion, where H spacc = 
span{|s)}™ =1 , % C oin = span{|0) , |1)}, and "H fus ion enu- 
merates all distinct measurement outcomes of topological 
charge when pairs of anyons are fused together [5] . 

One step of the walk is defined as a composition of 
two unitary operations, W — TU, where U = ^ ( \ \ ) 
acts in the coin space, and T is a conditional braiding 
operator which moves the walker left or right depending 
on the coin state: 

n 

T = I s - !) ( S M°) (0|<8>$.-i + |« + 1) (*|<8|1) (1|«4 , 

where b s = 6 S; i • • • 6 s , ms , b s = b s>mil ■ ■ ■ b s ,i, and b s = 
b s = 1 if m s = 0. The operators {b s ,i s }, acting on the 
fusion space % fusion) form a unitary representation of the 
r-strand braid group, r = 1 + X)l=i m sj which reflects 
the type of anyons we choose. To make T unitary, we as- 
sume periodic boundary conditions (|0) spaco = |«) spacc ) 
but will be concerned with walks satisfying n/2 < t, 
where t is the number of steps, so that winding around 
the surface is not an issue. 

Let the system's initial state be 1^(0)) = \sq) |cq) |$o)j 
where Sq = \n/2~\ is the initial position of the walker, 
c = denotes the initial state of the coin, and $o de- 
pends on the initial state of the anyons. After t iter- 
ations of the one step operator W, the state becomes 
|^(t)) = W l |^(0)) — a superposition over all coin histo- 
ries a € {0, 1}®*, weighted by appropriate phase factors. 
The reduced state of the spatial degree of freedom of the 
walker is 

pspacc (<) = tr coin tr fusion |*(0)(*(t)| 

= 53 trW S3' tr ^'|sS><S3'l , (2) 



where ss = s + J2k=i(^ a k — 1) ^ s tne walker's final po- 
sition corresponding to the coin history a = (ai, . . . , at); 
trUss- = j,(-l) z ^ s '\ with z(a,a f ) = Efc=i(ofcajk+i + 
a 'k a 'k+i)i i s a partial trace over the coin DOF; and 
yss' = Ba\$>o) (&o\Bt, acts in the fusion space. The 



braid word Bg can be constructed recursively from a 
given coin history a: 



Bjtk+i 



)bs s(k) -iB S (k) if afc + i — 
bs s(k) B S (k) if afc+i = 1 



(3) 



where = (oi, . . . , a&) is a truncation of a, s S (o) = So, 
and Bg(o) = 1. 

The spatial distribution of the walker after t steps is 
given by diagonal elements of the reduced density matrix, 



Prn(s,t) = (s\ Pspacc(i) \s) 



1 



(-ir^'hvyss, 



(a, a')- 



(4) 

where "(a, a') s" denotes the set of pairs of paths 
(a, a') satisfying at = a' t and = S3' = s. The subscript 
rfi indicates a fixed island occupation configuration. Vari- 
ance of the probability distribution p^{s, t) is defined in 
the usual way: cr?.(i) = Y, s Prn(s,t)s 2 - (Y, s P™(s,t)s) 2 . 

In the studies of transport phenomena in disordered 
environments, one is interested in quantities that result 
from averaging over all random background configura- 
tions. We shall assume that the island occupation num- 
bers m s are independent and identically distributed ran- 
dom variables with distribution W(m s ). The probabil- 
ity of occurrence of a configuration m is then simply 
Wm — n«=i W(m s ), and we denote configuration av- 
erage of a quantity Q A by ((Q)) = ^^A- Tnc 
average position distribution after a t-step walk is given 
by 



«p(M)» 



l 

2* 



(5) 



The topological quantity ((tr^gg/)), that depends on par- 
ticle statistics, governs the transport of an anyonic walker 
in random background. 

Abelian anyons:- For Abelian anyons the braid gener- 
ators {b s , is } are all equal to e 1 ^ . We assume the any- 
onic exchange angle to be 4> = ±j^, N € N. The fu- 
sion space is one-dimensional, % fusion — C, and we can 
choose |$o) arbitrarily. Upon introducing the linking 

, „ #(£„ and b B in Bg)-#(6+ and b+ in bL) 

numbers i s (a, a ) = ^ > 

that count the number of times the walker's trajec- 
tory (a, a') winds around an island s, trj^ss' reduces to 

For simplicity we consider uniform occupation distri- 
bution: W(m) = l/N for 1 < m < N. Then the av- 
erage position distribution of the walker after t steps, 
((p( ± ^)( S) t))) ; is given by (5) with 



«try£* J >> = nw* 



') mod N j 



(6) 



where Si j is the Kronecker symbol. Fig. 2 shows numer- 
ical results of ((p^ ±: ^\s,t))) for large time t and time 
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FIG. 2: Numerical results for localisation of Abelian anyons. 
The exchange statistics is <f> — ^ and the statistics is aver- 
aged over a random background of island occupations where 
the distribution in each is uniform over m 3 € {0, . . . , 7}. (a) 
Average variance as a function of time t for up to 600 time 
steps. The averages are taken over at least 500 charge con- 
figurations. (For clarity, only every 10th step is plotted.) (b) 
Average of the logarithm of probability distribution at time 
step t = 1000, taken over 10000 charge configurations. Red 
lines correspond to bounds for the localisation length given in 
Eq. (7). (Only the relevant region —700 < s — so < 700 and 
every 20th site are plotted.) In both a) and b), the error bars 
are given by standard deviation. 

behavior of the average variance 

((cr 2 (±f )(£)}). Further- 
more, we find that the results are essentially insensitive to 
the choice of N (except for the case N = 1 which corre- 
sponds to fermions, and which thus reduces to a usual 
quantum walk). The variance approaches a constant 
value and the asymptotic position distribution assumes a 
characteristic exponential shape, ((p^ ± w , )(s,t — > oo))) ~ 

l°-»0l 

e ? i°c , with the localization length £i oc = 1.44. This 
is a clear manifestation of the Anderson localization of 
Abelian anyons. 

Theoretical explanation of the above numerical results 
is based on the correspondence between our Abelian any- 
onic quantum walk in disordered background and the 
one-dimensional multiple scattering model presented in 
[10]. In the latter, scatterers are arranged in line with 
random distances between neighbours (see Appendix). 
Incoming light undergoes a series of scattering events and 
eventually localises due to the randomness in phases that 
individual trajectories accumulate during their passage 
between consecutive scatterers. The spatial randomness 
can be translated into a disorder in the population of 
the islands in the anyonic quantum walk model. Fur- 
thermore, reflection and transmission coefficients of the 
scatterers can be identified with the entries of the coin 
operator. The scatterer distances take discrete and fi- 
nite many values reflecting the discrete nature in the 
populations of the islands. Localisation of the walking 
Abelian anyon then follows, with the localisation length 
6oc ~ l/ln2 = 1.443. More precisely, in Appendix we 
derive the bound 

1 <r < 1 

In 2 + ln(l + 2^/2)2/^ - Voc - m 2 + hr(l - 2^/2)2/^ ' 

(7) 



where the anyonic statistical angle <fi = jt ■ This theoreti- 
cal result is in agreement with exact numerical treatment 
presented in Fig. 2. 

Ising non- Abelian anyons:- We consider the Ising 
model of non-Abelian anyons. There are three kinds 
of particles in this model: vacuum, fermion, and non- 
Abelian anyon. In analogy with [8], the initial fusion 
state |$o) describes the vacuum configuration of pairs of 
anyons with half the members braided to the right. Here, 
they randomly populate islands (except for the walker) 
to form a disordered background configuration m. The 
braid generators {b s ,i s } are now unitary matrices, the di- 
mension of which grows like d' m ', where |m| = 5Z™ =1 m s, 
and d = v2 is the quantum dimension of Ising anyons. 

The trace over the fusion degree of freedom can be 
related to the Kauffman bracket polynomial (Lss 1 ) [H] 
of a link Lggi, which arises from the Markov closure of 
the braid word Bl,Bg [6, 8, 12]. Moreover, (L$g>) can be 
expressed in terms of the Jones polynomial Vj,„, (q) [13]. 
Altogether, 

(8) 

where the writhe w(Lggi) = 2^2 s=1 m s £ s (a,d'). For 
Ising anyons, which correspond to spin-1/2 irreps of the 
quantum group SU (2)2 we have specifically q — i. In 
this case Vl ss , (i) can be further simplified in terms of 
simple topological characteristics of the link Lgg/ (see 
Appendix). Our final result for tr^gg' reads 

n 

was* = n u-i) i ~* ms n (-i)"™' r(sv '> , 

s=l Ks' <s" <n 

m s >0 

(9) 

where £ s = So,e 3 mod 2 and r(s',s") is the Milnor triple 
invariant of a three-component sublink of L^s' formed 
by strands corresponding to the walker and the islands 
s' and s" . (We do not explicitly indicate the (a, (in- 
dependence of the right hand side.) For the case of uni- 
form background configuration we recover the case stud- 
ied in [8]. 

To illustrate the role of disorder in a non-Abelian any- 
onic quantum walk, we choose uniform island occupa- 
tion probabilities. As (9) is for m s > 4-periodic in m s 
we take W(m) = 1/4 for 1 < m < 4 and W(m) = 
otherwise [14]. The configuration average of (9) can be 
partially carried through, as shown in Appendix. The 
average position distribution of the Ising quantum walk 
after t steps, ((p^ lsin6 '(s, t))), is given by (5) with average 
fusion space trace 

n 

((tr^£ ins) )) = Tm II SoMW) mod s , (10) 

s=l 

where Ta* - £ £* g{ o,i}» ni<.'< / '< n (-l) TBp,n - r( '''*" ) 
can be interpreted as an arithmetic mean of the quan- 
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FIG. 3: Numerical results for transport of Abelian and non- 
Abelian anyons. (a) Exact results for the discrete time quan- 
tum walk with anyons over an n = 46 sized lattice. Here 
is shown: variance of an Abelian anyon with | exchange 
statistics around a uniform background of islands singly occu- 
pied by Abelian anyons (Uniform Abelian=UA), variance of 
a Ising model anyon around a uniform background m s = lVs 
of Ising anyons (Uniform non-Abelian=UnA), variance of an 
Ising model anyon around a random background of islands 
where the distribution in each is uniform over m a £ {0, . . . , 4} 
Ising anyons and averaging is taken over 50 background con- 
figurations (Random non-Abelian=RnA) , and variance of the 
Abelian anyon averaged over at least 100 background config- 
urations of Abelian anyons with occupation probability uni- 
form over m s G {0, . . . , 7} (Random Abelian=RA). (b) Vari- 
ance of anyons in a Hubbard model on a ladder realising a 
continuous time anyonic quantum walk over an n = 100 sized 
lattice. Here space and time axes are scaled so that a contin- 
uous time classical diffusion would have diffusion coefficient 
providing a 2 (t) — 1. The asymptotic slope for the case RnA 
averaged in the same way as in (a) is 0.5426 and the slope 
remains positive and less than one within one sigma variance 
(error bars suppressed for clarity) . Numerics are obtained us- 
ing an approximate method employing real time evolution of 
an anyonic MPS with bond dimension equal to 100. 



tity (— l) ar ^( i aa') taken over all sublinks L* m , of a link 
Lggi. On comparing (10) to the Abelian expression (6) 
for N = 8, they are identical except for the prefactor 7ss' ■ 
Considering Eq. (6) as the coherent expression where the 
quantum interference of probability amplitudes causes lo- 
calisation, the Tss' coefficient can be viewed as a noise 
term which might or might not destroy the interference. 
In the Supplementary material, we argue that at short 
time scales, this term does not preserve memory and in- 
troduces temporal randomness. By results of Ref. [15], 
localisation does not occur in the presence of both spa- 
tial and temporal randomness, therefore we conjecture 
that non-Abelian anyons do not localise in the asymp- 
totic limit t — > oo. 

We have also performed numerical calculations for non- 
Abelian anyons using two methods. First, the probability 
distribution was calculated using Eqs. (4) and (9) up to 
23 time steps (see inset of Fig. 3), but it is hard to de- 
termine from these results whether localisation occurs or 
not. To obtain results for longer time scales, we used a 
Hubbard model [ ] on a ladder realising a continuous 
time anyonic quantum walk for an n = 100 sized lattice 



with open (reflecting) boundaries. These results were ob- 
tained by using the "Time-Evolving Block Decimation" 
(TEBD) [17] algorithm that is based on Matrix Product 
States (MPS). Fig. 3 shows that for non-Abelian anyons 
the variance grows linearly as a function of time, indicat- 
ing no signature of localisation. 

Conclusions:- A new characteristic of statistical be- 
haviour has been presented in terms of the localisation 
properties of anyons. We demonstrated that Abelian 
anyons propagating in an environment of other anyons 
of the same type in random positions, localise, regard- 
less of their mutual anyonic statistics. To prove this we 
demonstrated that localisation can occur not only in the 
presence of a continuous random variable, but also of dis- 
crete, where the localisation length is tightly determined. 
For the non-Abelian case we considered the Ising anyonic 
model. We showed that it fails to localise due to the en- 
tanglement of the walker with its environment that can 
be effectively described as dephasing of the walk leading 
to linear dispersion. This argument is supported by exact 
numerical results for a small number of walk steps and 
by an approximate numerical method that employs ma- 
trix product states. We expect delocalisation to occur 
even for more general SU{2)^ non-Abelian anyons due 
to the entanglement of the walker with the environment. 
Indeed, in [18] it has been demonstrated that the entan- 
glement of the walker with its uniform environment for 
any k > 1 is strong enough to cause the walk to decohere 
obtaining eventually a classical diffusive behaviour. 

Finally, our study provides a new paradigm of disorder 
that has the ability to override the localisation proper- 
ties of quasiparticles with anyonic statistics. In a real 
system with random potentials Abelian anyons that do 
not localise due to the potentials can be localised just by 
the presence of a random anyonic environment. On the 
other hand, non-Abelian anyons that would localise due 
to some random potential can be delocalised if a random 
background of anyons of the same type is present. 
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of Excellence in Engineered Quantum Systems (EQuS), 
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Multiple scattering model of an Abelian anyonic 
quantum walk 

Propagation of Abelian anyons and waves that undergo 
multiple scattering on a series of randomly distant scat- 
terers is governed by the same interference mechanism. 
We base correspondence between the two on physical in- 
tuition rather than mathematical rigor. 

The scattering model is described by Fig. 4. A 
monochromatic wave incident from the left scatters on 
a series of scatterers characterized by "from left / from 
right" reflection and transmission coefficients rj,tj/r'j,t'j. 
The distance between two successive scatterers j and j+1 
is random, such that the phase that the wave acquires 
when traveling between j and j + 1 is e l8 ' . Our result on 
localization of waves propagating through a random con- 
figuration of scatterers will be merely a discrete version 
of [1], section 2. 
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FIG. 4: In the multiple scattering model, the wave approaches 
(from the left) a series of n scatterers, and is transmitted with 
the amplitude The scatterers are arranged in line with 

random distances between neighbours. Hence, the phases e 1 j 
that the wave acquires during travelling from a scatterer j to 
j + 1 are also random. The complex quantities rj,tj /r'j,t'j 
are the reflection and transmission amplitudes for the wave 
impinging from the left / right. 

Denote by t\. n the block amplitude of transmission 
from the "left of scatterer 1" to the "right of scatterer 



n" ; and by r[ n the reflection amplitude from the block 
"1 to n" when approaching from the right. t\_ n can be 
expressed by the series 



k=0 

1 r n'i,n-l e 



(11) 



The corresponding transmission probability and its log- 
arithm are given by 



\tl,n-l\ 2 \t n \ 2 



ll-Vr^ 2 *"-*! 2 



(12) 



ln \t ltn \ 2 = ln |i 1 , n _ 1 | 2 +ln |i„| 2 -ln \l-r n r[ ^e™^ | 2 . 

(13) 

The reflection and transmission amplitudes t\ tn , r[ n are 
random variables that depend on the configuration of 
scatterers 1, . . . , n, i.e. on the angles Gi, . . . , 6 n -x- We 
assume that 9^s are identically distributed independent 
random variables with a uniform distribution over the 
discrete set {fjin \ m — 0, . . . , N — 1} (fj will be iden- 
tified with the anyonic exchange angle ip). We shall de- 
note by (((•••))) the statistical average over the angles 
8i, ... , 6 n -i, i.e. 



N-l 



N-l 



mi— m n —i=0 

Averaging of (13) leads to 
((l n |t 1>n | 2 » = ((ln|i 1 ,„_ 1 | 2 »+ln|t n | 2 



(14) 



- «ln|l- r n r' l n _ 



i26l„_i |2 



I 2 )} - (15) 



To proceed, we carry out the #„_i-average in the last 
term of (15), 



N-l 1 

E ^7 ln |l - r " r 'l.n-l e 



m n — 1=0 



I in 

N 



N-l 

n a 

m n _i=0 



•T m «-1 I 2 



, e N 



1 

N 



\n\l-(r n r') N \ 



N\2 



In the latter equality we used the fact that 

JV 



1-C 



N 



(16) 



(17) 



m— 1 



which can be proven by using the Newton's identities 
between elementary symmetric polynomials and power 



G 



The 6\, ■ ■ ■ , #„_2-average of (16) becomes trivial once 
we estimate (|r' 1 „_ 1 | < 1) 

Hl~(K\\r' M \) N ) 2 < mll-^r'^n 2 

< H^ + i\rn\\r' lin ^\) N )U8) 
Hl-\r n \ N f < lnll-^r'^n 2 

< Hl + \r n \ N ) 2 . (19) 

The upper and lower bounds of relation (15) read 

<(ln|ii,„| 2 )) < ((ln|t 1 ,„_ 1 | 2 ))+lnK„| 2 -lln(l-|r„| Ar ) 2 

(20) 

and 

<(ln|ti,„| 2 )> > ((ln|t 1 , Jl _ 1 | 2 ))+ln|i„| 2 -lln(l + |r„| Ar ) 2 

(21) 

respectively. When applied repeatedly, these recurrences 
yield (t 1A = ti) 



n 1 n 

«ln|i Ml | 2 ))< ^ln| £ ,| 2 --^ln(l 



\N\2 



J=2 



n i n 

<(lnK Ml | 2 )>> ^ln|^-| 2 -^ln(l + |r,r) 2 



J=2 



(22) 



where we have omitted the lower index of the averaging 
brackets ((. . .)). 

We shall now assume that tj — t,rj = r for all j. On 
the level of the Abelian anyonic quantum walk, this cor- 
responds to a spatially independent coin. Exponentiating 
(22) results in 



exp 
exp 



((Inl^l 2 ))^ (l-| r r ) i e -"H 1 -H N )*- 1 "l*l 2 ] , 

«m|t ln | 2 )}> (l + lrr)*^"^ 1 ^")*^ 1 !*! 2 ] . 

(23) 



Estimates of the localization length £; oc follow: 

1 , 1 

< Cloc < 



ln(l + Irl^)^ - ln|i| : 



ln(l — Irl^)^ — In |t| 2 
(24) 

For N -> oo we have £; oc -j^p- 

For the upper bound in (7) to make sense, — ln|t| + 
jj ln(l — H^) 2 has to be a positive number. This leads 
to the condition 

\t\ 2 < (1- \r\ N )^ , i.e. \t\ N + \r\ N < 1 . (25) 

Since |t| 2 + |r| 2 = 1 (with \t\, \r\ < 1), the latter is satisfied 
for N > 2. 

The case N — 1 corresponds to fermions which are 
known not to localize (their exchange statistics does 
not induce any interference effects). The marginal case 



N = 2 corresponds to semions (ip — ir/2), but we 
are unable to decide about their localization within this 
method. However, numerics for an analogous model in- 
volving continuous time hopping of semions on a ladder 
support localisation [ ]. 

To establish a connection between this scattering 
model and the Abelian anyonic quantum walk with the 
coin U = ( } \ ) , we define 



t 



1 



t' 



1 



1 



1 



(26) 



The localization length estimate for N = 8 (|-anyons) is 
now 



1.412 < &oc < 1.477 



(27) 



Let us stress that we investigated stationary state of a 
wave after infinitely many scattering events. This cor- 
responds to the infinite-time asymptotic behavior of the 
anyonic quantum walk. 



Simplification of the Jones polynomial Vl„, (i) 

The Jones polynomial for the links relevant to the 
quantum walk, can be related to a simpler arf invari- 
ant through [4] 



(28) 



s=l 
m s >0 



where £ s = 5o.i 3 mod 2- The product in the last expression 
is equal to 1 only if the link is proper, i.e. the sum of 
the pairwise linking numbers is even . Furthermore, when 
the link is totally proper, i.e. all pairs of components have 
an even linking number, and there is no self-linking, then 



arf{L sw ) 



s=l 



m s c 2 {s) + ^2 m s trn s "T(s',s") 



Ks' <s" <n 



where C2(s) is the cubic coefficient in the Conway poly- 
nomial of the two-component sublink of Lggr consist- 
ing of strands corresponding to the walker and island 
s; and t(s',s") is the Milnor triple point invariant of 
the three-component sublink of Lggi consisting of the 
walker's strand and the strands corresponding to islands 
s' and s" . Note there is no self- linking because the links 
represent world lines of the anyons which move forward 
in time only. In Ref. [(>] it was shown that for a uniformly 
filled background, for all the links that contribute to the 
diagonal elements of the spatial probability distribution, 
properness implies total properness. For non-uniform fill- 
ing the same is true for the simple reason that paths that 
had even linking between the walker and one island anyon 
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will then have multiple pairwise linking when the island 4 mod 8, then t(s', n) = for all r = 1, . . . , n — 1[18], 

has multiple occupancy. an d therefore (-1)^ (-l) E "'=i m =' r ( s '>") = -1. Hence, 

Inserting (28) and (29) into the expression for the fu- ((tty aa ,)) = whenever l n ^ mod 8. By the same 

sion space trace (8), we obtain reasoning, analogous holds for any island s = 1, . . . ,n. 

n Let us define i s = So i mo d 8- We have 

„?r>o ((tr^)) = Il4 E 4-"(-l)£i-'-"<" m =' m ^V') 

n s=l me{l,...4}» 

= (_!)£,'<." "V"V,r( s V') JJ 4(i)% m «(-l) m = c =( s ) . (35) 

The expression (— l) m a' m «" T ( s > s ) i s invariant under 
shifting m s i — > m s ' + 2 or m s >> — > m s >> + 2. There- 
fore, the sum over m £ {1, ... , 4}™ contains 2™ classes of 
2" equivalent configurations. Also, rrij — is equivalent 
with rrij = 2. We conclude 



s=l 
m s >0 



For islands s such that m s > and £ s 
we can use a result of [6] , and simplify 



1, i.e. 



(30) 



e z, 



02(a) = ^(£-1) = -- 
6 v s ; 2 3 



Hence 



- 1 



mod 2 <-s 



(31) 



n 



i 

2" 



iG{0,l}" l<s'<s"<n 

(36) 



trfe = n ^H)^ m * n i~i) m ^"< s '-^ . 

(32) 



s=l 
m 3 >0 



Ks' <s f, <n 



Configuration average of try^s' in Ising model 
anyonic quantum walk 

We assume uniform island occupation probabilities: 
W(m) — 1/4, 77i = 1, . . . , 4. The average trace over the 
fusion space for Ising model anyons is calculated as 



1 " ~ 



me(l,...,4)" s=l 



X 



(33) 



l<s'<s"<n 

We observe the following: 
1 

1 1 

s=l 



{(try*?)) 



1n-l 



e nc-o* 



mG{l,...,4}' 1 - 1 s=l 



l<s'<s"<n-l 
4 



. . . x 



i if (_i)^(_i)Er'=i m .' T (''' n ) = i 



otherwise 



(34) 



Hence, ((tr^ss')) = whenever £ n ^ mod 4 (i.e. when 

i 

(— i) 3 is not a real number). Furthermore, if £ n = 



Correlations in time 

To analyse the long-term behaviour of the non-Abelian 
case, we draw an analogy to decoherent quantum walks 
by interpreting the anyonic quantum walk as a quantum 
walk with an environment. The fusion space is a very spe- 
cial kind of environment, since it is non-local and highly 
non-Markovian, and such environments are not well stud- 
ied in the literature. Thus, anyonic environments open a 
new line of study on effects of environment and decoher- 
ence in quantum systems. 

Disordered quantum walks have been studied in the 
literature recently, and localization has been observed 
for many types of spatial disorder in the coin param- 
eter. Linden and Sharam [7] showed that periodically 
varying coin parameters can lead to both bounded and 
unbounded walks, depending on the period. Konno [8, 9] 
used path-counting methods to study quantum walks un- 
der quenched and annealed disorder, and showed that a 
coin defect at the origin causes partial localization of the 
wavepacket. Joye and Merkli [10] and Ahlbrecht et al 
[1 1] showed that dynamical Anderson localization occurs 
for spatially inhomogeneous coin, when the coin parame- 
ters are chosen randomly from continuous or certain dis- 
crete sets. If the coin parameters change in time, the 
walker spreads ballistically or diffusively, but no local- 
ization has been observed. Brun et al [12] studied peri- 
odically changing coins, and found that the walk is still 
ballistic unless the period becomes as large as the length 
of the walk, in which case it is diffusive. Shapira et al 
[1.3] studied quenched unitary noise numerically and ob- 
served that the walk is diffusive in the long term after an 
initial ballistic period. Similar conclusions were drawn 
in analytical treatment by Ahlbrecht et al [14] and Joye 
[ I- 
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As shown in Ref. [15], the quantum walk is diffusive 
also in the presence of both temporal and spatial disor- 
der, in other words localization does not occur in a spa- 
tially disordered system if temporal randomness is also 
present. The quantum coherent terms that are needed 
for localization in the spatially disordered case are thus 
destroyed by decoherence from temporal randomness. 
When a coin changes completely randomly, the identity 
of the next coin is not dependent on the previous coins at 
all, and the environment which induces the change has no 
memory. If the rate of change of the coin is larger than 
a single time step, there exists a Markov process charac- 
terized by some time- varying parameter "/(t). If there is 
no memory at all, these constants are identical for every 
t. In a Markov process, the future states of a system 
depend only on its present state and 7, not its history. 
In the following, we argue that the anyonic environment 
is effectively memoryless and is therefore described by a 
Markov process, which implies that the walker spreads 
diffusively. 

Equation (10) in the main text shows that the aver- 
age trace ((tr^ , l^, lng ' ) )) of Ising anyons is almost identical 
to that of Abelian anyons, except for the term involv- 
ing the Milnor triple invariant r. The Milnor triple in- 
variant gives the number of Borromean rings in a three- 
component link. We argue that the 7sa'-term induces 
temporal randomness in the walk, such that the en- 
vironment acts effectively in a Markovian way. The 
effect of this term is that it multiplies the contribu- 
tion from each path by the configuration average over 
(_1)£i< 8 '<s"<„™.'"V't(sV'). while this term preserves 

memory of the whole history of the particle's trajectory, 
we argue that at short time scales it fluctuates in a dis- 
ordered manner. The value of r changes when a new 
Borromean ring is formed, which requires at least 4 time 
steps. Also, the formation of a Borromean ring requires a 
very specific pattern in the particle's trajectory which is 
in no way periodic. In addition, because of the condition 
on the last coin outcomes at = a' t , there are new path 
patterns up to t — 1 time steps introduced on every time 
step which were not allowed for the previous time step, 
so these paths are not correlated to previous evolution at 
all. 

The time correlations of the r invariant can be tracked 
by defining the correlator 



C(t,t') = 



((-i)^(-ir^')-((-ir')((-i)w) 

i-((-iH 2 



where r t = Ei<s'<s"<ti t ( s 'j s ") is the sum of three- 
component invariants for all sublinks for a path up to 
t time steps and ( • )/gg/w> ao i s the expectation value 
over all paths leading to the initial site so (the subindcx 
has been suppressed above for clarity). The term in the 
denominator is a normalization factor, which is defined 
to be the value of the term in the numerator in the per- 




5 10 
Correlation time t' 



FIG. 5: Correlator C(t,t') as a function of t' with total num- 
ber of time steps t — 18. The line shows the best exponential 



fit C(i') = 1.7679 ( 



-0.57699 t' 



fectly correlated case t' = 0. For simplicity, the correlator 
is calculated only for the uniform filling (m s = 1 V.s) 
using the same method as described in [6]. The time 
correlations can be analysed by keeping the final time t 
fixed and calculating the correlator for increasing values 
of t' . The intermediate time value r t _ t / is calculated by 
erasing the t' last braid generators from the total braid 
word, such that braiding is switched off after t — t' time 
steps. For some braid words, the quantity r(s',s") is 
not well defined in our method, in which case we set 
t(s', s") = 0. The correlator for t = 18 is plotted in Fig. 
5, which shows that the correlations fall off exponentially 
after one step, indicating that effectively the r invariant 
maintains memory only for a short period of time, and 
the environment is Markovian at long time scales. Note 



that because of the condition a t 



the braiding at 



the last time step can always be trivially undone, there- 
fore the last two time steps are perfectly correlated. The 
correlator can also be calculated for other sites s 7^ sq. 
In these cases, the rapid falloff is observed for the central 
region ,s — < s < s + t/\/2, and outside this region 
the falloff becomes linear at the furthest edge sites. The 
edge behaviour is however irrelevant, as the behaviour 
of the quantum walk is determined by the central region 
only. 

The fundamental reason for the classical-like behaviour 
of the walker is the strong entanglement between the 
quantum walk states and the fusion states, and the highly 
mixing nature of the fusion space environment. The 
walker states stay entangled with the fusion states for 
long time periods, and recurrences where these states be- 
come uncoupled happen very rarely. 

We have also investigated the effect of temporally ran- 
dom phases introduced to the spatially random Abelian 
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walk. In this model the wave function is multiplied by a 
random -1 phase with some probability if the walker 
crosses a site belonging to the temporally fluctuating re- 
gion of sites. The calculations with different values of 
p_i and different sizes of the region up to 500 time steps 
showed that the behaviour becomes diffusive in all these 
cases. 

Numerical calculations of non-Abelian anyons 

We have calculated numerically the average variance 
((<7 2 ( Ising )(f))) over configurations m s € {0, ...,4} for 
the probability distribution given by Eqs. (4) and (9) up 
to 23 time steps (46 anyons in the lattice), with at least 
100 charge configurations. The error bars were obtained 
by the standard deviation of the variance of the spatial 
probability distribution. The average variance is approx- 
imately a straight line with slope 0.456 from 10 to 23 time 
steps, but the error bars of the variance overlap with the 
errorbars of the Abelian case, so we cannot distinguish 
between Abelian and non-Abelian anyons on this short 
time scale. Considering only occupations m s G {1, . . . , 4} 
(no vacuum charges), the variance is smaller, ie. the wave 
packet is diffusing slower (not shown). 

To obtain results for longer times, we turn to the 
continuous-time picture where the time evolution is gen- 
erated by a Hubbard-type Hamiltonian H and the propa- 
gator e - lHSt implements an infinitesimal time evolution. 
The total Hamiltonian is given by the sum of the shift 
and coin flip terms H — -ff s hift + -f^flip, where [16] 

#shift = jJ2( T r+APi + T+h s p 2 ) + h.c. 

s 

ff ap = 5>|2> a <l|+K*|l> a <2|) 

S 

withJeM,KeCandr s ± = (|l) s±1 s (l| + |2) s±1 ,<2|)<8 
^fusion are translation operators between sites s and s±l, 
b s are braid generators as defined in the main text for 
braiding the mobile anyone around anyons in the island 
between s and s + 1 (acting on fusion space only) and 
P c = J2 S l c ) s ( c l ® fusion are projectors to the coin states. 
Here \c) s corresponds to occupation of state c at site 
s, i.e. |c= 0) corresponds to no mobile anyon at site 
s, |c = l) s (|c = 2) s ) is a mobile anyon with coin state 
|0) (|1)) at site s. Here we consider only the case where 
coins are identical on every site. 



The above Hamiltonian is the generator of continuous- 
time evolution for total time T. Running the continuous- 
time walk for time T simulates the discrete-time quan- 
tum walk in a stroboscopic manner, such that the walker 
makes T/St steps of infinitesimal length St: e~ lHT — 
^ e -iH6tyr/5t _ j n first order of Suzuki- Trotter expan- 
sion, the propagator decomposes to e~ tHBMttSt e~ tHmpSt , 
similarly as the single step operator in the discrete-time 
quantum walk model. 

The TEBD algorithm was used to perform real time 
evolution of an initial state with the mobile anyon placed 
at the middle of the lattice and with couplings k = J = 1. 
Our implementation of the TEBD algorithm explicitly 
preserves anyonic charge [16] and also particle number 
[17] corresponding to the presence of a single walker. 
Thus, while the position of the walker and the fusion 
degrees of freedom of the anyons are entangled during 
the time evolution, the algorithm preserves distinction 
between them. This in turn allowed for specification of 
total anyonic charge and particle number, here we con- 
sidered the total vacuum sector with one walker. 
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